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Cristobalite is a common mineral in volcanic ash produced from dome-forming 
eruptions. Assessment of the respiratory hazard posed by volcanic ash requires 
understanding the nature of the cristobalite it contains. Volcanic cristobahte 
contains coupled substitutions of Al"*^ and Na^ for Si"*^; similar co-substitutions 
in synthetic cristobalite are known to modify the crystal structure, affecting the 
stability of the a and P forms and the observed transition between them. Here, 
for the first time, the dynamics and energy changes associated with the a-^ 
phase transition in volcanic cristobalite are investigated using X-ray powder 
diffraction with simultaneous in situ heating and differential scanning 
calorimetry. At ambient temperature, volcanic cristobalite exists in the a form 
and has a larger cell volume than synthetic a-cristobalite; as a result, its 
diffraction pattern sits between ICDD a- and /i-cristobalite library patterns, 
which could cause ambiguity in phase identification. On heating from ambient 
temperature, volcanic cristobalite exhibits a lower degree of thermal expansion 
than synthetic cristobalite, and it also has a lower a-fi transition temperature 
(~473 K) compared with synthetic cristobalite (upwards of 543 K); these 
observations are discussed in relation to the presence of Al^* and Na* defects. 
The transition shows a stable and reproducible hysteresis loop with a and P 
phases coexisting through the transition, suggesting that discrete crystals in the 
sample have different transition temperatures. 



1 . Introduction 
1.1. Background 

Cristobalite, a crystalline silica polymorph, forms in 
volcanic lava domes through vapour-phase mineralization and 
devitrification of groundmass (Baxter et ai, 1999; Horwell et 
al., 2013). Volcanic domes are mounds of highly viscous lava 
that pile up over the vent rather than moving away as a lava 
flow. In these domes, metastable cristobalite is found in 
preference to quartz as its formation is kineticaUy favoured 
(Brodie & Rutter, 2000), and because of the large activation 
energy required for the reconstructive phase transition from 
cristobalite to quartz. Cristobalite is stable at high tempera- 
tures and low pressures (1743-2001 K at 1 atm; 1 atm = 
101.325 kPa) in its cubic form, and undergoes a spontaneous 
displacive phase transition to tetragonal a-cristobalite upon 
cooling. The P to a transition is accompanied by a ~5% 
reduction in volume, which causes the crystals to crack on 
cooling (Carpenter et ai, 1998); the resulting 'fish-scale' 
texture is indicative of cristobalite that has undergone the 
displacive transition. 



As a result of the P to a transition, /S-cristobahte is not 
expected in nature under ambient conditions. However, 
several X-ray diffraction (XRD) studies of cristobahte in 
volcanic rocks have identified both the a and P forms. For 
example, Keith & Muffler (1978) analysed cristobalite in tuff 
deposits in the YeUowstone caldera, USA, and found both a- 
and j6-cristobalite. Keith et al. (1978) further identified 
/i-cristobahte in hydrothermally altered rhyolite flow deposits 
at Yellowstone. More recently, Reich et al. (2009) observed 
polycrystalline aggregates of cristobalite nanofibres in ash 
erupted from Chaiten volcano, Chile. These were identified as 
/i-cristobahte by diffraction pattern indexing of cubic lattice 
parameters a = b = c = 7.17 A, characterized by the 111, 220, 
331 and 422 reflections. 

Volcanic domes are inherently unstable and prone to 
collapse. On collapse, the cristobalite -bearing lava fragments 
to produce cristobalite-rich ash, prompting concern for 
exposed populations (Horwell & Baxter, 2006). Crystalhne 
silica is a known human carcinogen and long-term exposure 
may lead to silicosis, a chronic fibrotic lung disease (Beckett, 
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2000; lARC, 1997); consequently, much work has been carried 
out in recent years to establish whether volcanic cristobalite 
can cause disease (e.g. Jones & BeruBe, 2011; Wilson et al., 
2000; Horwell, Fenoglio et al, 2003; Forbes et al, 2003). 
Toxicological studies report low toxicity for cristobalite- 
bearing ash compared with a-cristobalite controls (e.g. CuUen 
et al., 2002) in short-duration experiments; however, the 
toxicity of y6-cristobalite is entirely unknown. 

The biological reactivity of inhaled crystalline siUca is 
governed by the origin and state (including whether freshly 
ground or aged) of the sample and its presence in a complex 
constitution, such as volcanic ash (Clouter et al., 2001; Fubini 
et al., 1995; Mossman & Glenn, 2013). As such, its pathogenic 
potential has been deemed a 'variable entity', whereby 
chemical and structural variability can modulate its toxicity 
(Donaldson & Borm, 1998; Napierska et al., 2010). The role of 
the crystallization environment and mineralogy on the toxicity 
of volcanic cristobalite has only recently been considered: 
Horwell et al. (2012) reported the substitution of up to 
2.4 mol.% AI2O3 for Si02 in cristobalite from Soufriere Hills 
volcano, Montserrat, and suggested that this could potentially 
reduce its toxicity since Al has been shown to quell crystalline 
silica toxicity (Duffin et al., 2001; Stone et al., 2004). In cris- 
tobalite, Al^"^ frequently occupies Si tetrahedral sites and Na"^ 
occupies the interstitial nonframework sites for charge 
balance (Perrotta et al., 1989; Thomas et al, 1994). Nonspecific 
structural modifications have also been implicated in the 
bioreactivity of crystalline silica (Guthrie, 1997), although 
never considered for volcanic cristobaUte. 

In this work, we exploit the a-fi phase transition to deter- 
mine the effect of the substitutions described by Horwell et al. 
(2012) on the crystal structure, and compare the behaviour of 
volcanic cristobahte with the known structural and thermo- 
dynamic character of synthetic cristobalite. This is a necessary 
step towards understanding the formation and preservation of 
metastable cristobalite in volcanic systems. Additionally, since 
the physicochemical properties that govern the biological 
response to crystalline silica are directly related to the origin 
of the sample (Fenoglio et al., 2000), this investigation furthers 
our understanding of the relationship between intrinsic 
mineralogy and the surprisingly limited in vitro toxicity of 
cristobalite-bearing volcanic ash. 

1.2. The a—fi phase transition in synthetic cristobalite 

The transition between the a and /3 forms of cristobalite is 
first order; hence, it is expected to be instantaneous and to 
correspond to a discontinuous jump in entropy (Schmahl et al. , 
1992). However, many studies of synthetic cristobahte samples 
have shown that the transition occurs over a temperature 
interval of ~20-40 K (Leadbetter & Wright, 1976; Mosesman 
& Pitzer, 1941; Sosman, 1965). This transition range exists 
because individual crystallites within a cristobalite-bearing 
sample may have considerably different transition tempera- 
tures depending on their size, structure and composition. 
Furthermore, a proportion of grains may persist metastably 
after they have passed their transition temperature. Conse- 



quently, at a given temperature within the transitional range, a 
and y6 forms coexist both stably and metastably (Leadbetter & 
Wright, 1976; Swainson et al., 2003; Schmahl et al, 1992; 
Spearing et al., 1992). The a-f) transition has been shown to 
depend on temperature, but not on time, demonstrating that 
the finite width of the transition range is not simply a kinetic 
effect (Schmahl, 1993; Leadbetter & Wright, 1976). A marked 
hysteresis is also consistent with the first-order character. 

The a-/S transition temperature in synthetic cristobalite is 
highly variable, ranging from 443 to 543 K (Frondel, 1962). 
This variabihty predominantly results from structural substi- 
tutions, along with the associated defects, vacancies and 
strains. Within this range, lower a-jS inversion temperatures 
are associated with structural imperfections; for example, 
when a mineralizer [e.g. Al2(S04)3 and Na2C03] is present 
during crystallization (Chao & Lu, 2002&; Stevens et al., 1997). 
It is known that cation substitutions in cristobalite cause a 
change in the lattice dimensions, which results in a shift in the 
prominent 101 a-cristobalite peak position observed by XRD. 
Chao & Lu (2002fl) have shown that this displacement is 
proportional to the level of (AI2O3 + Na20) co-doping during 
cristobalite synthesis; they find that dioi increases from 
4.040 A at 2.91 mol.% Al203-Na20 to 4.058 A at 4.65 mol.%. 
These levels of doping are of the order of those seen for 
volcanic cristobalite (up to 2.4 mol.% AI2O3). The a-/3 tran- 
sition may be suppressed entirely when substitution is exten- 
sive (e.g. Chao & Lu, 2002i>; Perrotta et al, 1989). However, in 
the Si02-Na20-Al203 system, stabilization of yS-cristobalite at 
room temperature was not observed below 6.29 mol.% AI2O3 
in synthetic samples (Chao & Lu, 2002fl). The fitoa transition 
in the closely analogous quartz system can also be suppressed 
by incorporation of other elements into the Si02 structure: Xu 
et al. (2001) found that the high-temperature p form of quartz 
can be preserved down to room temperature through the 
substitution of Al^^ and Li* for Si"**. In this case, suppression is 
due to a reduction in the transition temperature to below 
ambient temperature. The presence of foreign ions in the 
interstices and altered bond angles from Al coordination 
presumably inhibit the contraction of the structure. 

The a-yff transition temperature and transition interval may 
also depend on the crystallization and annealing history of the 
sample, since these influence crystallite properties. Lower 
transition temperatures are associated with lower crystal- 
lization temperatures; for example, Alcala et al. (1996) observe 
a transition temperature of 523 K for Si02 annealed at 1673 K 
versus 503 K when aimealed at 1373 K. By carefully control- 
ling the conditions of synthesis, it is possible to create cristo- 
balite samples in which the a-fi transition interval is very 
narrow. Thompson & Wennemer (1979) report a sharp peak in 
heat capacity at 535 K for a highly crystalline synthetic cris- 
tobalite (prepared at 1773 K for 20 h), indicating that a narrow 
transition range may reflect a homogeneous crystal popula- 
tion. Conversely, multiple phases of disordered growth and 
incorporation of impurity ions could increase the breadth of 
the transition interval (Stevens et al., 1997). 

The a-fi transition temperature can also be a function of 
crystal grain size, which is dependent on duration of forma- 
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tion. In synthetic samples, increased annealing times lead to 
larger average grain sizes, although grain growth occurs 
heterogeneously (Lee & Lee, 2000). Analogously, in the 
quartz system, the thermodynamics of the a-fi transition differ 
depending on the origin and size fraction of the sample; 
changes in the primary 101 quartz peak intensity that often 
accompany the a-fi quartz inversion are a function of particle 
size (Wahl et al, 1961). In general, large particles favour the 
spontaneous formation of low-temperature (a) phases, 
whereas small particles stabilize the high-temperature (/i) 
phase (Lee & Lee, 2000), and phase transitions can become 
sensitive to grain size in nanoscale crystals owing to the 
relationship between atomic positions and macroscopic strain 
(Pertsev & Salje, 2000). Rios et al. (2001) show this for the 
transition from y6 to a quartz, but the change from normal 
properties to 'nano' properties, when the bulk properties 
become dominated by surface effects, occurs below "-'50 nm 
(McKnighte? a/., 2008). 

2. Materials and methods 

2.1 . Volcanic ash and cristobalite reference samples 

Bulk ash sample MBA5/6/99 from Soufriere Hills volcano, 
Montserrat, West Indies, was analysed (Montserrat Volcano 
Observatory reference number MV01163). It was collected 
during a dome-collapse event on 5 June 1999 and has been 
used extensively in ash characterization studies {e.g. Horwell, 
Fenoglio et al, 2003; Horwell, Sparks et al, 2003). The sample 
contains ~11 wt% cristobalite (Horwell et al, 2014). 

Two cristobalite standards were analysed as controls: a 
pure-phase synthetic sample (denoted DKSmith) and an 
impure synthetic sample (denoted Hemenway). The samples 
were supplied by Deane K. Smith (Pennsylvania State 
University) and David Hemenway (University of Vermont). 

2.2. Concentration of crystalline silica component 

The sub-4 (im 'respirable' fraction of MBA5/6/99 (denoted 
MRA5/6/99) was also analysed. It was separated from the bulk 
sample using a British Rema Minisplit Aerosol Classification 
System (Process Plant Engineers, Sheffield, UK) at 
16 000 r min^'^ to concentrate the cristobalite, which is found 
in greater abundance in the finer component (Horwell, Sparks 
et al, 2003). 

Crystalline silica was further concentrated from the respir- 
able ash sample MRA5/6/99 by the Talvitie method (Talvitie, 
1951), in which all minerals other than crystalline silica are 
dissolved by boiling the sample in phosphoric acid at 293 K for 
8 min, and residues removed by washing with fluoroboric and 
hydrochloric acids. XRD revealed that the concentrated 
crystalhne silica fraction was composed of 97 wt% cristobalite 
and 3 wt% quartz, although a trace of plagioclase feldspar (the 
dominant mineral phase in the bulk sample) remained, with 
perhaps also a trace of tridymite (Horwell, Sparks et al, 2003). 
The residual plagioclase is likely to be due to the intimate 
association between cristobalite and other phases in individual 
ash particles (Horwell et al, 2012). 



2.3. Characterization of cristobalite and phase transitions 

Cristobalite in all five samples was identified using X-ray 
diffraction with a position-sensitive detector (XRD-PSD) at 
the Natural History Museum, London. This instrument 
consists of a Nonius X-ray diffraction system with an INEL 
curved 120° 2Q PSD with the beam and sample in static 
reflection geometry. We use the cristobalite peak of highest 
intensity (near 22° 26) for identification in our unconcentrated 
volcanic ash samples because lower-intensity cristobalite 
peaks are often masked by the presence of other minerals: 
plagioclase, in particular. XRD patterns were compared with 
the PDF-2 database cards from the ICDD (International 
Centre for Diffraction Data, http://www.icdd.com) for a- and 
/S-cristobahte for identification purposes. 

To date, the ICDD database contains 14 theoretical and 
empirical powder diffraction patterns for cristobalite with 
sufficient information for comparison. These library patterns 
were overlaid using WinX''°'^ (STOE, Darmstadt, Germany) 
in order to construct a database range for the 101 a- and 111 /S- 
cristobahte peaks (see Fig. 1). The ICDD PDF-2 cards used 
for ambient a-cristobahte were 00-039-1425, 01-071-0785, 01- 
075-0923, 01-076-0935, 01-077-1316, 01-082-1232, 01-089-3434 
and 01-089-3606. Two further PDF-2 cards for a-cristobalite 
were included: 01-082-0512, which is a crystal structure 
calculation for a-cristobalite with distorted ions, and 01-076- 
0941, which is a pattern collected at 503 K. The PDF-2 cards 
used for ;6-cristobalite were 01-085-0621, 01-089-3435, 01-089- 
3607 and 01-076-0934. 

The a-/j cristobalite phase transition in samples MBA5/6/ 
99, MRA5/6/99 and Hemenway was investigated by XRD 
using a GeniX Cu high-flux X-ray system equipped with a 
Xenocs FOX two-dimensional 10_30P mirror to generate 
high-brightness Cu Ka radiation (Natural History Museum, 
London). Samples were heated in situ, using an Anton Paar 
heating stage (HTK-10) with a platinum strip, from room 
temperature to 573 K with a stepwise heating profile at a rate 
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Figure 1 

Range of (101)„ (red) and (111)^ (blue) cristobalite peak positions 
defined by ICDD library patterns. Vertical lines correspond to the peak 
positions of the 14 ICDD patterns under consideration. Dashed lines 
indicate patterns *01-082-0512 (structure with distorted ions) and **01- 
076-0941 (503 K). The diffraction profile of sample MBA5/6/99 is shown 
to demonstrate where volcanic cristobalite typically sits within these 
ranges. 



j. Appl. Cryst. (2014). 47, 1205-1215 



David E. Damby et a/. • Phase transition in volcanic cristobalite 1 207 



research papers 



of 10 K min^\ The high-brightness X-ray source allowed 
rapid data collection of isothermal diffraction patterns, each of 
only 1 min duration, at 2 K intervals. Thus, a sequential 0-120° 
W snap-shot visualization of the changes in diffraction peak 
positions throughout the heating cycle was recorded by the 
PSD. The point at which the (lOl)o. peak began decreasing in 
intensity marked the onset temperature of the a-yS transition, 
and the point at which the (IOI)q, peak became indis- 
tinguishable from background defined the end temperature. In 
order to monitor the effect of prolonged heating on the a-jS 
transition and (lOl)o, d-spacing shift, additional heating 
experiments were performed, where samples were held at 
various mid-transition temperatures. 

The position of the (101)^ cristobalite peak is affected by 
the overlap of a minor plagioclase peak, which appears as a 
shoulder on the cristobalite peak at higher 2Q angle; however, 
the plagioclase peak does not overlap with the (111)^ cristo- 
balite peak. Thermal expansion of plagioclase with composi- 
tions typical of Soufriere Hills (An>5o; Couch et al, 2003) is 
relatively limited over the temperature variation in the 
present study, and most of the volume expansion is accom- 
modated for anisotropically in the direction of the (100) plane 
normal (Tribaudino et al, 2010); therefore, the plagioclase 
peak near 22.0° 29 was considered to be unaffected by heating 
to 573 K. Consequently, it could be isolated from patterns 
collected at temperatures above the a-P transition and 
subtracted from the lower-temperature patterns, allowing the 
(101)„ cristobaHte and plagioclase peaks to be deconvolved. 
Samples were left undisturbed during the entire heating 
sequence so that patterns for different temperatures could be 
directly compared. Since the plagioclase peak position 
remained fixed through consecutive runs, this is a good indi- 
cation that no change in sample-stage height/sample surface 
occurred during the thermal heating experiments; therefore, 
no internal 29 calibration standard was required and accurate 
peak positions were achieved using an external Si standard 
alone. Peak positions for the primary (101)„ and (111)^ cris- 
tobalite reflections are discussed according to the 29 position 
at half-width half-maximum. 

Differential scaiming calorimetry (DSC) was used to 
provide a complementary determination of the a-fi transition 
temperature for sample MBA5/6/99 observed by XRD. These 
results provide a transition-temperature range unobstructed 
by the plagioclase overlap as well as associated transition 
enthalpies for volcanic cristobalite. Samples were heated from 
423 to 573 K (encompassing the entirety of the a-^ transition 
range) on a PerkinElmer Pyris 1 differential scanning calori- 
meter in the Department of Chemistry, Durham University. 
The thermal stability of samples was assessed by thermo- 
gravimetric analysis (TGA) on a PerkinElmer PjtIs 1 ther- 
mogravimetric analyser in the Department of Chemistry, 
Durham University. 

2.4. Rapid quenching of volcanic )?-cristobalite 

Samples of volcanic ash were heated through the cristoba- 
lite transition and cooled at different rates to investigate 



Table 1 

Summary of d spacings and transition parameters for samples 
investigated. 

Lattice data are in angstrom and temperatures are in kelvin. £/(101)„ and 
are tlie d spacings for primary a- and /i-cristobalite peaks, respectively. 
Arf„ is tlie thermal expansion of rf(101)„ for tlie a pliase prior to tlie onset of 
the a-fi transition. Adc^f) is the change in the d spacing during the transition 
following thermal expansion (Ad„). T,, is the transition temperature range for 
the a-/S inversion detected by XRD and DSC. 

rf(ioi)„ d(ni)/, r„- r„ 

Sample Type 298 K 573 K Ad^ Ad^^p (XRD) (DSC) 

MBA5/6/99 Ash 4.064 4.124 0.014 0.046 443-483 448-508 

MRA5/6/99 Ash 4.067 4.133 0.015 0.051 448-493 - 

Hemenway Standard 4.055 4.108 0.018 0.035 513-553 - 

DKSmith Standard 4.038 - _ _ _ 

Talvitie Isolate 4.057 - _ _ _ 



whether rapid cooling could preserve j8-cristobaUte. MBA5/6/ 
99 was pressed into 10 mm pellets at 9.8 MPa, heated from 
ambient temperature to 573 K at 5 K min~^, and held for 1 h 
before being either rapidly quenched in liquid nitrogen or left 
to cool at room temperature. Pellets were gently ground and 
run immediately by XRD on a Bruker AXS D8 ADVANCE 
with DAVINCI design in 29 reflection mode using Cu radia- 
tion and an Ni filter (Department of Chemistry, Durham 
University). Data were interrogated using the DIFFRAC. 
EVA software (Bruker AXS Inc., Madison, WA, USA). A 
control pellet was pressed (but not heated) and analysed by 
XRD under the same conditions. 

2.5. Imaging and elemental composition of volcanic cristo- 
balite 

Cross sections through cristobahte particles (MBA5/6/99) 
were prepared for analysis by mounting a sub-sample in resin 
(Araldite resin AY103 and hardener HY956, 5:1 ratio). The 
resin block was coated with 20 nm of carbon. Imaging and 
elemental analysis of cristobalite crystals were undertaken in a 
Hitachi SU-70 field emission gun scanning electron micro- 
scope equipped with an Oxford Instruments INCAx-act (LN2- 
free analytical silicon drift detector) energy-dispersive X-ray 
spectroscopy (EDS) system at the G. J. Russell Microscopy 
Facility, Durham University. Imaging was used to confirm the 
presence of the characteristic 'fish-scale' cracking in cristo- 
balite and EDS was used to qualitatively confirm the presence 
of cation substitutions observed by Horwell et al. (2012). 

3. Results 

3.1 . Identification of volcanic cristobalite using library 
patterns 

A comprehensive comparison of ICDD library patterns for 
cristobalite revealed two distinct ranges for the a and forms 
(Fig. 1). The primary (101)„ peak at ambient conditions is 
between dioi = 3.982 and 4.047 A (21.9 to 22.3° 29). When 
considering the non-ambient (503 K) and distorted a-cristo- 
balite patterns, the (101)„ peak range extends to dioi = 4.081 A 
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Degrees 29 

Figure 2 

XRD patterns of the Hemenway cristobalite standard (left) and ash sample MRA5/6/99 (right) at four distinct temperatures: (a) room temperature as 
a-cristobalite; (b) prior to the a-P transition following thermal expansion (Arf„); (c) midway through the transition; and (d) following the transition to fl- 
cristobalite. Patterns were collected for intervals of 1 min per given temperature and are displayed between 21.0 and 22.5° 26 to enlarge the primary 
(101)„ and (111)/! peak region. 



(21.8° 29). The primary (111)^ peak is between = 4.082 
and 4.137 A (21.5 to 21.8° 20). 

The d spacings observed at ambient temperature in this 
study for volcanic cristobalite were 4.064 (MBA5/6/99) and 
4.067 A (MRA5/6/99), and 4.057 A after the Talvitie treat- 
ment (Table 1), which fit within the range for a-cristobalite 
identified from the ICDD patterns but are larger than the 
range compiled for a-cristobalite under ambient conditions. 
No evidence of /J-cristobalite was observed in samples of 
volcanic ash at ambient temperature. The two synthetic cris- 
tobalite standards also fit within the range for a-cristobalite 
identified from the ICDD patterns: dun = 4.055 A (Hemen- 
way) and dioi = 4.038 A (DKSmith). 

3.2. Characterization of the displacive transition 

Superimposed thermal XRD patterns at temperatures 
below, across and above the a-P transition from continuous 
runs for both volcanic cristobalite and the Hemenway 
synthetic sample show a clear change in the position of the 
primary cristobalite peak with increasing temperature, 
marking the a-P transition (Fig. 2). The transition occurs over 
a temperature range of around 40 K, taking place between 443 
and 483 K in sample MBA5/6/99, 448 and 493 K in sample 
MRA5/6/99, and 513 and 553 K in the Hemenway sample. 
DSC data show a consistent, but slightly broader, temperature 
range of 448-508 K for the transition in volcanic sample 
MBA5/6/99. 

Thermal XRD data show that the cristobahte transition 
comprises a gradual migration of the (101)„ peak position 
(A(i„) towards lower 26 angles, followed by a definitive change 
in d spacing at the a-fi transition {Ada_^) (Fig. 2). For volcanic 
cristobalite, Ada_p ranges from 0.046 to 0.051 A. This change 
in primary peak position follows a A^/„ of 0.08° 26 (0.014 A) 
from room temperature to the onset of the transition. These 
phenomena were also observed in the Hemenway sample, but 
with a smaller Ad^-p of 0.035 A following a larger migration of 
Ad„ of 0.018 A. 



DSC data for MBA5/6/99 show two clear endothermic 
reactions during the first heating cycle, one at 403 K and 
another at 475 K (Fig. 3). Only one exothermic reaction took 
place on cooling, near 456 K, and only the higher-temperature 
transition recurred on reheating, indicating that the peak at 
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Figure 3 

(fl) The ff-jS transition of cristobalite detected by DSC (sample MBA5/6/ 
99) with an initial heating cycle, a cooling cycle and a second heating 
cycle. The peaks at ~475 K on heating and ~456 K on cooling mark the 
phase transition, (b) Heat capacity of the anomalous peak centred on 
~403 K and corresponding mass loss measured by TGA. Scanning rate is 
10 K mm-\ 
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403 K is a nonreversible event (discussed in §4.5). As no otlier 
minerals identified in volcanic ash go through a transition in 
this temperature range, the mirrored reactions at 475 and 
456 K are recording the transition from a- to /S-cristobalite on 
heating and the reverse on cooling. The temperatures at which 
the phase transition occurs on heating and cooling do not 
correspond, however. The onset temperature of the transition 
during one cycle corresponds to the maximum transition 
temperature of the other cycle; i.e. the a to yS transition on 
heating commences at the peak temperature of the cooling 
transition (449 K onset versus 456 K cooling cycle maximum), 
and vice versa for the /i to a transition (475 K for both tran- 
sition onset and heating cycle maximum) (Fig. 3). The tran- 
sition shows an increased integrated heat flow for a 
consecutive structural rearrangement over two scans, from 
9.001 to 10.030 J g"^ for MBA5/6/99 (Fig. 3). The susceptibility 
of cristobalite to thermal hysteresis explains the difference in 
the transition onset temperature upon heating and cooling. 

The a-P phase transition proceeds as a function of 
temperature rather than time, as demonstrated when the 
transition was paused and the sample held isothermally mid- 
transition. An isothermal step maintained for 1 h upon heating 
through the transition range failed to progress the phase 
transition further; however, an additional increase in 
temperature following the isothermal step immediately 
resumed the inversion. When the temperature was held 
constant, the d spacing remained the same as that initially 
measured during the regular heating cycle. Additionally, the 
same proportion of a- and /J-cristobalite was observed to exist 
at discrete temperatures in repeated experiments. 

Cooling rate had no effect on the final crystal structure after 
first heating MBA5/6/99 to produce the p phase. No differ- 
ences were observed in the cristobalite diffraction patterns for 
the non-heat-treated sample, the rapidly quenched sample and 
the sample slowly cooled to ambient conditions. 

No change in mass was observed by TGA through the a-fi 
transition range. However, a mass loss of ~0.12 wt% was 
observed from 383 to 413 K (Fig. 3). 

The peak positions for the Talvitie-treated separate differ 
from the corresponding non-treated sample MRA5/6/99 {d^ = 
4.057 and 4.064 A, respectively). This probably results from 
the boihng of the sample within the a-fi transition tempera- 
ture range in the Talvitie method. The sample was not 
subjected to thermal analysis since the treatment appears to 
have affected the sample's thermochemistry, and further 
examination is, therefore, outside of the aims of this study. As 
a result, it is unknown whether intraparticle phases in volcanic 
ash identified by Horwell et al. (2012) impose a constraint on 
the cristobahte transition. However, data are included in 
Table 1 to show the effects of the concentration process on the 
crystal lattice. 

3.3. Observation of cristobalite in ash particles 

The distinctive 'fish-scale' texture was observed in cross 
sections through volcanic cristobalite (Fig. 4). Occasionally, 
individual crystals within ash grains did not display the 'fish- 
scale' texture, but these were always less than ~8 |j,m in 



diameter and probably represent either a minimum crystal size 
for which the change in volume results in shattering or indi- 
vidual 'fish-scales' of cristobalite crystals that fragmented 
along the cracks during dome collapse. A critical grain size of 
>5 |j,m has been observed for which the thermally induced 
'fish-scale' cracking can exist (Lee & Lee, 2000), which is 
appropriate for the present observations. Scanning electron 
microscopy (SEM)-EDS confirmed the presence of alumi- 
nium and sodium in cristobalite previously observed in sample 
MBA5/6/99 by Horwell, Sparks et al. (2003). The abundance of 
substitutions was not quantified in the present study, as 
previous analyses of cristobalite from Soufriere Hills have 
shown a range of 0.3-2.4 mol.% AI2O3 within distinct samples 
(HorweU et al., 2012). 

4. Discussion 

The presence of cristobalite with 'fish-scale' cracking in 
volcanic ash (Fig. 4) provides unequivocal evidence that it 
crystalhzed in the dome as /S-cristobahte but subsequently 
transformed to the a form following eruption and cooling. 
Cristobalite with diagnostic cracking has also been observed in 
dome rock from Soufriere Hills (Horwell et at, 2013; Plail et 
al., 2014), Merapi (Damby et al., 2013), Mount St Helens 
(Blundy & Cashman, 2001; Pallister et al., 2008) and Santia- 
guito (Damby, 2012; Scott et al., 2012), and sectioned ash from 
Sakurajima (Hillman et al., 2012). In the following discussion 
we present evidence from experiments by thermal XRD 
(Fig. 2) and DSC (Fig. 3) on volcanic ash from Soufriere Hills 
to explain the dynamics of the total (excepting nanofibres, see 
§4.6) transition of fi- to a-cristobahte during cooling. 

4.1. Structure of volcanic cristobalite 

The interplanar spacings calculated by XRD for volcanic 
a-cristobahte (d^oi = 4.064-4.067 A) (Table 1) are larger than 




Figure 4 

Backscattered electron image of a cristobalite-rich particle from ash 
sample MBA5/6/99 in polished section. Patches of cristobalite marked as 
'C are two 'fish-scales' characteristic of the 'fish-scale' cracking texture 
that occurs as cristobalite transforms from the /i phase to the a phase on 
cooling. Lighter-grey patches are feldspar and glass. Image taken at 
8.0 kV and working distance of 15.0 mm. 
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for the Hemenway sample (rfioi = 4.055 A), for the DKSmith 
sample (rfioi = 4.038 A), the range defined by ICDD patterns 
for a-cristobalite at ambient conditions and reported values 
for synthetic a-cristobalite with few impurity ions (e.g. diQ^ = 
4.04 A; Chao & Lu, 2002fl); they trend towards the (101)„ peak 
positions of the non-ambient and distorted ICDD a-cristo- 
balite patterns. The larger d spacing observed in natural 
a-cristobalite probably results from the structural substitu- 
tions observed by Horwell et al. (2012). However, at elevated 
temperatures, the interplanar spacings of volcanic j6-cristo- 
balite {dm = 4.097-4.124 A) correspond well to those in the 
Hemenway sample at elevated temperature (dm = 4.108 A), 
as well as literature values for synthetic /3-cristobahte (dui = 
4.11 A) (Butler & Dyson, 1997) and the presently defined 
ICDD /6-cristobalite range. 

The peak widths for volcanic cristobalite are wider than 
those for the two synthetic cristobalite standards used, indi- 
cating a larger degree of disorder in the volcanic crystals and 
variable short-range interplanar spacings. A homogenous 
distribution of substituted cations would be expected to shift 
the 29 diffraction peak profile whilst maintaining half-peak 
width to height ratios; however, Damby (2012) reports a 
heterogeneous distribution of substituted ions throughout 
individual cristobalite crystals. Such chemical variation will 
result in non-uniform d spacings and lead to a broader XRD 
peak profile, as is observed, and supports the conclusion of 
structurally incorporated cations in volcanic cristobalite. 
Poorly ordered cristobalite results in an XRD pattern that has 
a higher backgroimd and decreased sharpness and intensity of 
reflections compared to a well ordered sample (Alcala et al. , 
1996); therefore, we suggest that structural disorder needs to 
be carefully considered when quantifying mineral phases by 
peak height [e.g. the IAS method of Le Blond et al. (2009)]. 

Caution should also be taken when determining the 
presence of a versus P polymorphs in volcanic ash on the basis 
of primary peak matching with the ICDD database. Such 
matching could potentially prove to be ambiguous (Fig. 1), 
which we attribute to three factors: (1) the increased lattice 
spacing of (Al,Na)-bearing cristobalite in volcanic ash 
samples, which broadens and shifts the (lOl)^, diffraction peak 
towards the position where the (111),3 peak occurs; (2) the 
ubiquitous presence in cristobalite-bearing volcanic ash of an 
overlapping plagioclase peak near 22° 26, which masks the 
true position of the cristobalite peak; and (3) the range in the 
ICDD peak positions reported from different laboratories. 
These factors can lead to the primary cristobalite peak in 
volcanic ash often appearing closer to certain library patterns 
for /S-cristobalite rather than a-cristobalite. Additionally, it is 
possible that user-based differences in sample preparation, 
instrument deployment and instrument calibration methods 
can result in different peak positions that are of variable 
accuracy. 

Notwithstanding these issues, the ranges for a and P forms 
defined by the ICDD library patterns may help with identifi- 
cation by eliminating the need to deal with the aforemen- 
tioned considerations. For example, in the study by Keith & 
Muffler (1978), which identified both a- and ;S-cristobalite in 



tuff deposits, the (101)„ peak was quite sharp and ranged from 
21.8 to 22.0° 26icu (fi^ioi = 4.07^.04 A), whereas the (111);, 
peak was usually rounded and ranged from 21.6 to 21.8° 20cu 
(dill = 4.11-4.07 A). These ranges are in good agreement with 
the present study. 

4.2. Characterization of the a—fi cristobalite transition 

The cristobalite inversion involves a transition enthalpy, 
which is recorded by the endo- and exothermic reactions 
observed by DSC (Fig. 3). No other major phase in these 
samples identified by XRD undergoes a first-order transition 
between 443 and 513 K, leaving cristobalite as the only 
discrete phase inversion in this range. The diffraction intensity 
and lattice-spacing data also confirm a structural inversion of 
cristobalite in Soufriere Hills volcanic ash (Ad„_^ = 0.046- 
0.051 A). The transition-temperature ranges observed by DSC 
(448-508 K) and XRD (443-483 K) are in close agreement 
and probably differ because of the different instrument 
heating/temperature measurement techniques and the manner 
in which transition temperatures were determined from the 
raw data. 

Comparing XRD patterns from prior to and after a heating 
and coohng cycle demonstrates that there is no change in 
cristobalite peak intensity or in the proportions of other 
phases present. Therefore, from the characteristics of the 
(IOI)q, peak, volcanic cristobalite reverts back to the original 
disordered isomorph upon cooling. This reversion of disor- 
dered a-cristobalite to the same defect structure following a 
transition was also observed by Butler & Dyson (1997) and 
probably results from the inability of foreign ions and 
vacancies that stabilize the disordered form to be substantially 
varied on the timescale of the experiments. Further, repeat 
analyses on sub-samples show that the same proportions of a- 
and j8-cristobalite coexist at a specific temperature. Therefore, 
the increase in the integrated heat flow from subsequent runs 
of the same sample (Fig. 3) does not necessarily imply an 
increase in cristobalite abundance or inherent structural 
order, but instead is likely to result from the irreversible 
nature of the transition hysteresis for cristobalite, where 
cycling through the a-/i transition is accompanied by an 
increase in entropy (Richet et al., 1982). 

The mixed-phase nature of volcanic ash makes it difficult to 
relate the heat of transition for the a-fi inversion to that 
observed for pure-phase synthetic samples. However, in 
scaling up the enthalpy associated with the transition in 
sample MBA5/6/99, which has ~11 wt% cristobalite on the 
basis of an XRD quantification by HorweU et al. (2014), the 
heat flow for a 'pure-phase' sample of volcanic cristobalite 
would be approximately 5400 J mol~^. This result is in poor 
agreement with values from previous studies on synthetic 
samples, of 2320 and 1300 J moP^ (Thompson & Wennemer, 
1979; Mosesman & Pitzer, 1941), especially since less energy is 
required to transform a disordered cristobalite sample (with a 
larger unit cell) compared with a better ordered sample 
(Butler & Dyson, 1997). We note, however, that the quanti- 
fication by Horwell et al. (2014) has an imcertainty of ±3 wt%. 
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which is problematic given that a difference of 1 wt%, from 11 
to 10 wt%, would result in ~500 J mol^^ difference in heat 
flow. This effect scales with the proportion of cristobalite and, 
therefore, the error in calculated heat flow will increase with 
increasing cristobalite abundance. Further, as a mixed-phase 
dust, sub-samples of volcanic ash are expected to contain 
variable quantities of mineral phases; therefore, it is not 
possible to confirm the abundance of cristobalite in the sub- 
sample used for the DSC measurements. As such, more 
accurate quantification of cristobahte in heterogeneous dusts 
would be required to constrain the effects of structural defects 
on the transition energy. 

4.3. Controls on the a—fi transition temperature 

The a-^ transition temperature for Soufriere Hills ash 
occurs over a broader temperature range than is generally 
reported for synthetic cristobalite samples: ~40 K for ash 
versus ~20^0 K for synthetic samples (Leadbetter & Wright, 
1976; Mosesman & Pitzer, 1941; Sosman, 1965). We propose 
that the breadth of the transition reflects heterogeneity 
amongst cristobalite crystal grains, resulting from the presence 
of substituted ions and associated structural defects. Cristo- 
balite crystallization can commence within hours or days of a 
magma package entering the dome environment (Williamson 
etal, 2010; Horwell et al, 2014), and crystals will be subjected 
to variable temperatures, pressures and interaction with 
hydrothermal fluids throughout their growth and residence. 
As a result of this dynamic thermochemical environment, 
volcanic cristobalite probably undergoes multiple phases of 
disordered growth and ion substitution. As discussed in §1.2, 
these factors can increase the transition breadth compared 
with a homogenous cristobalite crystal population (Thompson 
& Wennemer, 1979; Stevens et al, 1997). 

The transition temperature for the cristobalite in our ash 
samples (443-483 K) falls at the low end of the range reported 
for synthetic cristobalite samples (§1.2) and is much lower than 
that in our Hemenway sample (513-553 K). We find that 
volcanic cristobalite exhibits a lower degree of thermal 
expansion (Ado,) than synthetic cristobalite, but expansion for 
both samples begins at similar temperatures: ~403 K for 
volcanic cristobalite and 393 K for the Hemenway sample. 
Consequently, we hypothesize that the lower a-^ inversion 
temperature observed in volcanic cristobalite may result from 
a smaller A.da associated with thermal expansion prior to the 
transition to ;6-cristobalite. This phenomenon would be 
facilitated by substituted and interstitial cations in volcanic 
cristobalite. O atoms in tetrahedra surrounding an AP^ 
substitution wiU be under-bonded and any such O atoms also 
adjacent to an interstitial Na* (satisfying local charge balance) 
are likely to occupy a slightly larger thermal ellipsoid and be 
dynamically (or statically) displaced from their preferred 
tetrahedral geometry within such a defect region. If the 
thermal energy input was transferred with higher probability 
to the O atoms local to impurity (Al^^, Na^) defect regions 
relative to other O atoms in the framework, then this may 
result in an increased thermal ellipsoid vibration for the O 



atoms incorporated at defects. In such a situation, the 
absorption of thermal energy primarily by increased oxygen 
vibrations in defect regions together with a more limited 
overall lattice expansion may be a sufficient response to cope 
with an incremental increase in thermal energy. Conversely, 
the lattice of a-cristobalite with few substitutions/defects may 
behave more uniformly and expand evenly to a greater overall 
degree in response to applied heat. In this way, local expansion 
and oxygen displacement in defect regions may trigger the a-P 
transition by lowering the activation energy. 

Cristobalite in the ash sample is likely to have formed 
through both vapour-phase mineralization and groundmass 
devitrification (Baxter et al, 1999; Horwell et al, 2012); a 
potential for these two mechanisms to produce crystals with 
different Al abundances has been reported (Damby, 2012). 
The above discussion would therefore suggest a potential 
difference in the a-yS transition temperature for cristobalite of 
vapour-phase and devitrification origin. However, we cannot 
test this hypothesis here because it was not possible to obtain 
separate samples for vapour-phase and devitrification cristo- 
balite. 

In synthetic cristobalite samples, lower transition tempera- 
tures are associated with lower crystallization temperatures 
(§1.2). The volcanic dome environment is typically below 
~1073 K (Barclay et al, 1998; Murphy et al, 2000), which is 
more than 300 K cooler than the minimum crystallization 
temperatures used to synthesize cristobalite; for example, a 
minimum temperature of 1373 K was used in experiments by 
Alcala et al. (1996) and Stevens et al (1997). By assuming a 
crystalhzation temperature of less than 1273 K in a volcanic 
setting, and extrapolating the observed effects in synthetic 
samples annealed between 1273 and 1673 K to cooler 
temperatures, we hypothesize that the lower inversion 
temperature seen for volcanic cristobalite may partially reflect 
the crystallization temperature in the dome. 

4.4. Coexistence of a and fi phases 

It is apparent from the continuous DSC endo- and 
exotherms (Fig. 3) and the gradual emergence of the (111)^3 
peak and disappearance of the (101)„ peak during thermal 
XRD (Fig. 2) that a and j6 forms coexist during the phase 
transition. This observation agrees with studies on synthetic 
cristobalite (e.g. Leadbetter & Wright, 1976; Swainson et al, 
2003; Schmahl et al, 1992; Spearing et al, 1992), whereby a 
smooth change in the intensity of these modes is observed as 
the proportion of cristobalite existing as either a or P varies 
with temperature (§1.2). The coexistence interval is larger in 
volcanic samples than has been observed for synthetic 
samples; two potentially complementary factors could 
contribute to the observed coexistence interval (§1.2): cation 
substitution and crystallite size. 

The range of cation substitutions observed in volcanic 
cristobalite crystals, from 0.30 to 2.38 mol.% AI2O3 (Horwell 
et al, 2012), may support the existence of different transition 
temperatures amongst grains, as the amount of ;S-cristobaUte 
retained metastably in experimental studies is known to 
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depend on the level of additives (Chao & Lu, 2002ft). Domains 
with a higher degree of disorder/abundance of additive ions 
are likely to transform at lower temperatures than purer 
regions. Since the temperature and duration of heating in 
these experiments are insufficient to allow significant diffusion 
of additives, it might be expected that the same transformation 
behaviour, in terms of coexisting aip proportions and the 
ranges of transformation temperature observed on cooling 
and re-heating, will be replicated exactly if the impurity 
regions are those where the transformation {via oxygen 
displacement) is initiated. 

The broad temperature range of the a-P transition could 
also result from polydisperse crystal sizes in our samples, 
where individual crystallites transform at different tempera- 
tures (§1.2). A heterogeneous cristobalite grain size distribu- 
tion has been quantitatively observed for volcanic dome 
samples, where larger crystals can form in domes that have 
longer residence times (Damby, 2012). A grain size control 
may be manifest in the small difference between XRD and 
DSC transition ranges as only the XRD samples were ground 
in a mortar and pestle prior to analysis. 

4.5. Anomalous DSC peak 

The endothermic reaction in the initial DSC heating cycle at 
403 K is a nonreversible event as no corresponding event took 
place upon cooling (Fig. 3). Most likely, this event is a dehy- 
dration reaction, resulting either from surface hydration and 
open porosity or from the crystallization of an amorphous/ 
hydrous phase such as opal-C or opal-CT. These two 
phenomena are indistinguishable in the current data as 
powdered opal requires minimal or no diffusion of 'free' water 
to the surface and, therefore, water is lost at a temperature 
below that of a corresponding lump sample (Smallwood et al. 
2008). Opal-C can be indistinguishable from cristobalite by 
XRD since the characteristic 101 peak around 21° 29 (4.00- 
4.10 A) of these materials coincides {e.g. Onal et al, 2007), 
although opal-C usually produces a broad diffraction peak 
compared with cristobahte. Considering a AH of vaporization 
at 403 K of 39.1 kJ mor\ this corresponds to 0.24 wt% of 
water in the sample, which is larger than the 0.12 wt% 
observed by TGA. This theoretical mass loss could correspond 
to as much as 5 wt% opal-C in the sample. The difference 
between the theoretical water content and experimental TGA 
data could be explained by a variation in mineral abundance 
in sub-samples from a single sample: we frequently observe 
variability of a few weight per cent for other phases in volcanic 
ash. This amount of opal-C (or opal-CT) would be near 
background level in the thermal XRD patterns as 60 s 
isothermal data collection snapshots resulted in peaks with 
intensities of only ~1000 counts. 

Alternatively, the discontinuity observed during the first 
heating cycle may reflect a tridymite inversion, if this phase is 
present at low abundance (not detectable by XRD) in the 
cristobalite. k transitions in tridymite are observed near 288 K 
in synthetic tridymite-cristobalite mixed phases (Thompson & 
Wennemer, 1979; Stevens et al, 1997). However, as all tridy- 



mite phase transitions are reversible (Graetsch & Florke, 
1991), a discontinuous inversion from tridymite to cristobalite 
would be required as no corresponding structural rearrange- 
ments are seen on further heating and coohng cycles. The 
change in latent heat observed in the first and second heating 
cycles (9.001 J g^^ for the first heating cycle versus 10.03 J g~^ 
for the second) could be interpreted as evidence of the 
transformation of tridymite to cristobalite, with incorporation 
of its associated latent heat into the subsequent cristobalite 
inversion. 

4.6. Preservation of ;S-cristobalite and implications for 
exposure 

From textural studies (Horwell et al, 2013) and the heating 
experiments contained herein, cristobalite crystallizes and 
persists in the volcanic dome in the P form when above the a- 
P transition temperature. However, during cooling, the 
charge-coupled substitutions in volcanic cristobalite are 
insufficient to suppress the transition to a-cristobalite as is 
possible with chemically stabilized cristobalite in the Na20- 
Al203-Si02 system (e.g. Chao & Lu, 20026; Perrotta et al, 
1989). As discussed in §1.2, a maximum of 2.4 mol.% AI2O3 
was observed in volcanic cristobalite, whereas stabilization of 
yS-cristobalite in synthetic samples required a minimum of 
6.29 mol.% AI2O3. Furthermore, the observed crystals are 
much too large for nanoscale effects (§1.2) to suppress the p to 
a transition for the bulk of cristobalite in volcanic ash. 
However, we propose that the preservation of /6-cristobalite 
nanofibres observed by transmission electron microscopy in 
ash from the Chaiten volcano, Chile (Reich et al, 2009), is an 
effect of grain size, since the diameter of the fibres (~20- 
50 nm) is below the length scale (~50 nm) at which nanoscale 
phase-transition effects have been observed (McKnight et al, 
2008). 

These observations, combined with our failure to preserve 
/S-cristobalite in the ash through rapid quenching experiments, 
imply that cristobalite in volcanic ash is overwhelmingly Al- 
bearing a-cristobalite following an eruption; hence, this is the 
phase of relevance to human exposure. From the present data, 
we conclude that a-cristobahte standards are appropriate to 
assess the toxicity of volcanic ash from a polymorphic 
perspective, although perhaps not with regards to structure (or 
composition, as shown by Horwell et al, 2012). However, if a 
small proportion of nearly amorphous hydrous silica is 
responsible for the anomalous DSC endotherm, this may have 
ramifications for toxicity as the bioreactivity will differ 
between cristobalite and other arrangements of Si02 {e.g. 
Pavan et al, 2013). 

5. Conclusion 

This is the first structural and thermodynamic investigation of 
volcanic cristobalite and exploits the a-P phase transition to 
determine its characteristics. Volcanic cristobalite conmionly 
exhibits substitutions, mainly Al^^ and Na^ for Si'*^ (Horwell et 
al, 2012), which serve to increase the size of the unit cell and 



/. Appl. Cryst. (2014). 47, 1205-1215 



David E. Damby et a/. • Phase transition in volcanic cristobalite 1213 



research papers 



to lower and broaden the a-^ transition temperature relative 
to more pure Si02 samples. Both thermal XRD and DSC 
demonstrate that a-cristobalite exists in volcanic ash, indi- 
cating that these substitutions are insufficient to stabilize 
j6-cristobalite at ambient temperature, as has been achieved 
with higher levels of Al and Na co-substitutions in chemically 
stabilized synthetic cristobalite. However, the observed 
structural alteration can lead to confusion when indexing 
against ICDD hbrary patterns owing to a shift in the primary 
(lOl)a peak to lower angles. The data presented here may help 
to clear up complexities shown by inverted samples and the 
broad range in the ICDD database; to facilitate this, we have 
defined distinct d-spacing ranges for a- (21.8-22.5° 20) and 
j6-cristobahte (21.5-21.7° 20) that could be used in the future 
for identification. 

Upon heating, volcanic cristobalite exhibits a smaller 
degree of thermal expansion (Ad^) of the (101)„ planar 
spacing prior to the a-f) transition compared with synthetic 
samples. Since the onset of this thermal expansion occurs at 
similar temperatures for both volcanic and synthetic cristo- 
balite, we suggest that the smaller Ad^ in volcanic cristobalite 
may be a direct result of the degree of substituted and inter- 
stitial ions, which also effectively cause the a-/3 transition to 
occur at lower temperatures. Owing to peak overlaps in 
mixed-phase volcanic ash, three-dimensional expansion could 
not be considered. No differences in ;S-cristobalite cell 
dimensions were observed compared with standard cristoba- 
lite samples or accepted values from studies of synthetic 
cristobaUte. 

Since crystalline silica from different environments has 
been shown to be a variable respiratory hazard in industrial 
settings, constraining the volcanic cristobalite hazard is 
imperative for hazard response during an eruptive episode. 
The verification of a-cristobaUte in volcanic ash means that 
toxicity studies which compare volcanic cristobalite with 
a-cristobalite standards are valid. However, natural cristoba- 
lite will always show some cation substitutions, and the hazard 
posed by exposure to respirable volcanic cristobalite should be 
considered accordingly. 
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